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Multi-band near-infrared images of twelve BL Lacertae objects were obtained with the 2.5m telescope at the Las 
Campanas Observatory in order to determine the properties of their underlying host galaxies. Resolved emission 
was clearly detected in eight of the lowest redshift targets (up to z ^--^0.3), and was modeled with a de Vaucouleurs 
r'/"^ surface brightness law. We find that the morphologies match the elliptical galaxy profiles well, and that the 
BL Lac objects reside in large and luminous, but otherwise normal hosts - consistent with previous studies done 
predominantly at optical wavelengths. The median absolute K-band magnitude of the galaxies in this study is 
-26.2, the average half-light radius is 4.2 ±2.3 kpc, and their average integrated R-K color is 2.7 ±0.3 mag. 
These are well within the range of values measured previously in the H-band by Kotilainen et al. and Scarpa 
et al. in a comparable number of targets. Taking their data together with our results, we find a best-fit K-band 
Kormendy relation of /i^. = 4.3 logio (r^/kpc) +14.2 mag arcsec"^, virtually identical to that obtained for normal 
ellipticals. Finally, the near-infrared colors determined for five galaxies (average J-K = 0.8 ±0.3 mag) are the 
first such measurements for BL Lac hosts, and match those expected from old stellar populations at the BL Lac 
redshifts. 



Subject headings: BL Lacertae objects — galaxies: elliptical and lenticular, cD 
parameters — galaxies; structure 



-galaxies: fundamental 



L INTRODUCTION 

The study of the host galaxies of BL Lacertae objects in re- 
cent years has contributed much to our understanding of these 
enigmatic objects. According to unified schemes for active 
galactic nuclei (AGN), BL Lacs are low-luminosity Fanaroff & 
Riley (1974) type I (FRI) radio galaxies whose jets are aimed 
very close to our line of sight (e.g., Urry & Padovani 1995). 
This explains their observed extreme properties primarily by 
relativistic beaming in the jet. The host galaxies and environ- 
ments of BL Lacs should therefore be independent of orienta- 
tion effects and should be similar to those of the hosts of FRIs. 
Indeed, observational data confirm that BL Lac objects are pre- 
dominantly found at the center of giant, but otherwise normal, 
elliptical galaxies (e.g., Wurtz, Stocke, & Yee 1996; Urry et 
al. 2000; Scarpa et al. 2000a). Also, the available integrated 
colors and color profiles in a number of these galaxies reveal 
stellar populations typical of normal ellipticals at the appropri- 
ate redshifts (Kotilainen, Falomo, & Scarpa 1998b; Urry et al. 
1999; Scarpa et al. 2000b), consistent with the bulk of the star 
formation having occurred at the redshifts of at least a few. 

Here, we report on moderately deep near-infrared (NIR) ob- 
servations of twelve BL Lacs drawn from the larger sample of 
110 objects observed (mostly in R-band) in the Hubble Space 
Telescope {HST) WFPC2 survey of BL Lac objects (Scarpa et 
al. 2000a; Urry et al. 2000; Falomo et al. 2000). The survey tar- 
gets were drawn from seven complete radio, optical, and X-ray 
flux limited samples. A subset of the sample was observed with 
NICMOS at H-band (Scarpa et al. 2000b), and there are two ob- 
jects in common with the present observations (MS 2143+070 
and H 2356-309). Our relatively long ground-based NIR ob- 
servations had the advantage over the NICMOS data that we 



were able to follow the surface brightness profiles of each ob- 
ject out to larger radii. NIR filters give us greater sensitivity 
to older, and therefore redder stellar populations compared to 
optical filters. The number of targets in this study is compa- 
rable to the H-band samples studied by Kotilainen, Falomo, & 
Scarpa (1998b) and Scarpa et al. (2000b), and we utilized their 
results for comparison. Following much of the previous litera- 
ture on this subject, we use Hq = 50 km s"' Mpc"' and = 
throughout this paper. 

2. OBSERVATIONS AND DATA REDUCTION 

Twelve BL Lac objects were observed with the du Pont 2.5m 
f/7.5 telescope at the Las Campanas Observatory on four nights 
from 28 September to 1 October 1998 using the 256x256 sq. 
pixel NIR camera (Persson et al. 1992) at medium resolution, 
with a scale of 0".348/pixel. The seeing during our observa- 
tions was typically below 1" FWHM (Table 1). We used the Jj 
(1.24 pm), H (1.65 /im), and K^ (2.16 pm) filters, where and 
K^ are variants of the respective standard bands designed to re- 
duce the background signal at these wavebands (Persson et al. 
1998), and H is the H-band filter of the CIT/CTIO photomet- 
ric system. Differences between the Jj- and K^-bands from the 
standard J and K are negligible and we refer to only the stan- 
dard bands throughout the rest of this paper We observed two 
BL Lacs in all three NIR bands, three were observed in both J 
and K, and the remaining seven were imaged in K only. The 
observations are summarized in Table 1 . 

We used the NOAO Image Reduction and Analysis Facility 
(IRAF^) package to calibrate the data. Raw exposures were 
dark subtracted and flat fielded. Flats for each filter were made 
by median combining several night sky exposures - usually re- 
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ferred to as "super sky flats". Bad pixels were identified via 
a mask image made from a ratio of two real sky flats with 
different exposure times using the task CCDMASK. These bad 
pixels were substituted with interpolated values using the task 

FIXPIX. 

Due to the high and widely varying sky brightness at NIR 
wavelengths, we dithered our exposures over nine positions in 
a square grid on the CCD to construct the sky frames, usually 
acquiring several short exposures at each position in order to 
prevent saturation. Each dithered image was sky subtracted, 
shifted relative to the center image, and then median-combined. 
In most cases, a series of these dithered exposures of the target 
were taken consecutively in the same night and were combined. 
Observations of standard stars taken from Persson et al. (1998) 
were used to set the photometry. We found an inconsistency 
in the photometry for one star (SS9175) used during the last 
night of our observing run. The chart for this object is sus- 
pected to have been printed incorrectly (E. Persson, 2000, pri- 
vate communication), therefore we simply did not use it in our 
calibration. The RMS in the zero-points derived from the three 
other standards within each night was smaller than 0.1 mag in 
all three bands. We found also that the standards were stable to 
better than 0. 1 mag over the course of Ihc four night observing 
run; we therefore estimate conservatively that our la error in 
photometry is 0.1 mag. 

In order to maximize the signal from the underlying 
hosts in our analysis, we extracted azimuthally averaged 
one-dimensional surface brightness profiles from the final 
BL Lac images using an interactive script implemented in 
SUPERMONGO. The profiles were then fitted with a point spread 
function (PSF) plus de Vaucouleurs (1948) r'/"* elliptical galaxy 
model convolved with the PSF, using a minimization tech- 
nique as outlined in Scarpa et al. (2000a). In all cases (with the 
possible exception of H 1914-194) the hosts were known to be 
elliptical (Scarpa et al. 2000a), so no attempts were made to fit 
the resolved emission with exponential disk profiles. 

Stars within our ^ 1.5' x 1.5' field of view were often too 
dim to extract a PSF with accurate wings. Since only the cen- 
tral part of the PSF is affected by the seeing, we combined the 
innermost (~ 2") part of the field star profile (or an average 
when multiple stars were in the field) in each BL Lac image 
with the outer wing extracted from a standard star observed in 
the same band on the same night. We found that this technique 
provided us with rehable descriptions of the PSF out to ~ 8". 

Final K-band images of the twelve objects observed are 
shown in Figure 1; the available 1- and H-band images display 
very similar morphology to the K-band images so are not pre- 
sented. Azimuthally averaged K-band profiles for the targets 
extracted from the images in Figure 1 are presented in Figure 2, 
and the best fit model parameters for all bands are reported in 
Table 2. In several cases, we observed the same object in the 
same band on two separate nights during our observing run. 
We analyzed these data independently and found that the de 
Vaucouleurs plus PSF fits on the data from the different nights 
agreed well within our la hmits. 

3. RESULTS 

The host galaxies of the eight BL Lacs at z < 0.308 were 
clearly resolved in our images, as they were also in the HST 
optical survey (Scarpa et al. 2000a; Urry et al. 2000). The 
three objects at larger (z > 0.5) redshifts were unresolved 
(PKS 0138-097, MS 0419+194 and PKS 2240-260) so only 



some brief notes are presented for them in section 3.3. In the 
HST survey data, the host of HI 9 14-1 94 was resolved but was 
consistent with both elliptical and disk morphologies. It is also 
resolved (though marginally) in the present data, but there is 
no robust determination of its redshift in the literature and fur- 
ther discussion of this object is hmited to section 3.3. The de- 
rived properties of several host galaxies studied previously in 
the H-band by Scarpa et al. (2000b, MS 2143+070, H 2356- 
309) and Kotilainen, Falomo, & Scarpa (1998b, PKS 0521- 
365, PKS 2005-489, MS 2143+070) are in agreement with our 
data (section 3.3). 

3.1. Morphologies and Luminosities of the Host Galaxies 

The resolved hosts are large, with average half-light radius 
(re) 4.2 ± 2.3 kpc. Although the uncertainties on the half- 
light radii are large (with the exceptions being the nearby 
bright galaxies in PKS 0521-365 and PKS 2005-489), the de- 
rived sizes are consistent with those derived from similar NIR 
ground-based observations by Kotilainen, Falomo, & Scarpa 
(1998b) where they found < >= 4.8 ± 1.9 kpc (omitting 
their large value of = 32.5 kpc for PKS 2254+074). On 
average, our results are also consistent within the errors with 
< r<, >= 10 ± 5 kpc found in the H-band NICMOS study by 
Scarpa et al. (2000b), although the ground-based data seem to 
prefer smaller values. A direct comparison of the derived 
for the two objects (MS 2143+070, H 2356-309) common to 
both the ground and space-based studies show that the individ- 
ual values are well in agreement (see section 3.3). 

Recently, Barth, Ho, & Sargent (2003) noted that the half- 
light radii derived from the HST R-band survey (Scarpa et al. 
2000a; Urry et al. 2000) were systematically smaller than those 
obtained through ground-based optical studies in their low red- 
shift sample of eleven BL Lacs. This difference is most pro- 
nounced in the largest galaxies where they suggest a combina- 
tion of short exposure time and relatively small field of view of 
the survey data may have led to an underestimation of re values 
from the HST. We do not find such large discrepancies in our 
comparison of the NIR data above - our comparison actually 
suggests on average that the situation may be slightly reversed 
from what is found in the comparison of optical data, so this 
issue remains unresolved. 

BL Lac host galaxies have been found to be quite luminous 
in optical studies (Falomo 1996; Wurtz, Stocke, & Yee 1996; 
Urry et al. 2000), with luminosities commensurate with those 
of the brightest cluster galaxies (Taylor et al. 1996; Wurtz, 
Stocke, & Yee 1996). Our median rest-frame galaxy NIR 
magnitudes, Mk = -26.2 mag (with a range of -26.9 to - 
25.6) and Mj = -25.4 mag, are about one magnitude brighter 
than the characteristic magnitude for nearby non-active ellip- 
ticals (M^ = -25.0 mag; Mobasher, Sharpies, & Ellis 1993). 
These observations are consistent with the H-band values found 
by Kotilainen, Falomo, & Scarpa (1998b) and Scarpa et al. 
(2000b), and overall, confirm results from the previous opti- 
cally based studies that the hosts of BL Lac objects are large 
and luminous, but otherwise normal, elliptical galaxies. 

Despite the fact that we were plagued with large uncer- 
tainties in our measurements of the effective radii, we com- 
bined our K-band /i^ - data (Table 2) with the H-band data 
from Kotilainen, Falomo, & Scarpa (1998b) and Scarpa et al. 
(2000b) to construct a K-band Kormendy relation for BL Lac 
host galaxies. The H-band data were converted to K assuming 
H-K = 0.22 (Recillas-Cruz et al. 1990). This Unear relation is 
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Fig. 1. — K-band contour images of the targets observed after smoothing with a gaussian filter with a= 1 pixel. In each case, a ~ 35" X 35" section of the field 
is shown, centered on the BL Lac object. Minimum contour levels (in mag arcsec ) are indicated on each image and successive isophotes are spaced by one mag. 
North is up on each image. For the objects which we obtained images in two separate nights, we only show the deeper image obtained with better seeing. 
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Fig. 2. — K-band azimuthally averaged radial surface brightness profiles of the sources (shown with error bars). These profiles were extracted from the images 
presented in Figure 1. Also shown are the scaled PSF (dotted line), the de Vaucouleurs r'/"* galaxy profile convolved with the PSF (dashed line), and their sum (solid 
line). Wherever only the solid line is shown, the host galaxy was unresolved. 
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a projection of the fundamental plane relating the half-light ra- 
dius to the surface brightness at that radius (Djorgovski & Davis 
1987; Kormendy & Djorgovski 1989). As a result, we obtained 
a best linear fit of /Xe = 4.3 logio (r^/kpc) +14.2 mag arcsec"-^ 
to the data (see Figure 3). This relation is virtually identical 
to that derived from K-band measurements of a sample of 59 
normal elliptical galaxies (/Xg = 4.3 logio (''e/kpc) +14.3 mag 
arcsec"^, after converting to our adopted cosmology) by Pahre, 
Djorgovski, & de CarvaUio (1995), and is consistent with the 
relation found previously using the H-band data only (Scarpa et 
al. 2000b). Thus, BL Lac hosts and normal elliptical galaxies 
appear to be dynamically similar and occupy the same region of 
at least this projection of the fundamental plane. Lastly, we fit a 
similar /ig-r^. trend in both slope (~4.4) and intercept ('--^13. 5) 
with a fairly large scatter using the resolved and marginally re- 
solved data for a similarly small sample of flat spectrum radio 
quasars (Kotilainen, Falomo, & Scarpa 1998a). This appears at 
least qualitatively similar to our derived relation for BL Lacs, 
indicating perhaps that the luminosity of the central source does 
not strongly influence the properties of their surrounding host 
galaxies. 
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Fig. 3. — K-band lie -re Kormendy relation for BL Lac host galaxies us- 
ing our K-band results (stars) and the H-band data from Kotilainen, Falomo, 
& Scarpa (1998b, triangles) and Scarpa et al. (2000b, circles) which were 
converted to K assuming & H-K = 0.22 color typical of normal eUipticals 
(Recillas-Cruz et al. 1990). The solid line is a linear least-square fit to all of 
the data (/ie = 4.3 logio (J'e/kpc) +14.2 mag arcsec"^) which is virtually identi- 
cal to that derived for a sample of 59 inactive elliptical galaxies in the K-band 
(dashed line: He = 4.3 logio (''e/kpc) +14.3 mag arcsec"^, after converting to 
our adopted cosmology) by Pahre, Djorgovski, & de Carvalho (1995), and is 
consistent with a H-band only relation of iXe =3.8 logio (''e/kpc) +14.8 mag 
arcsec"^ obtained previously by Scarpa et al. (2000b). 

3.2. Host Galaxy Colors 

We calculated integrated rest-frame 7?-/^ colors for the 8 of 
the 9 galaxies resolved in our observations (neglecting H 1914- 
194, which has no measured redshift) utilizing R-band values 
derived from the uniformly analyzed HST survey data (Scarpa 
et al. 2000a; Urry et al. 2000). The results are reported in Ta- 
ble 3. Additionally, we have information in two NIR filters for 

* Using <V-R>= 0.6 mag (Gregg 1989) aad<V-K >= 3.3 mag (Schombert 
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five of these galaxies, so we were also able to calculate J-K 
colors. To our knowledge, these are the first NIR colors pub- 
lished for any BL Lac host galaxy. We obtained average R-K 
and J-K values of 2.7 ± 0.3 and 0.8 ± 0.3 mag, respectively. 
These colors closely match those found in normal non-active 
elliptical galaxies: < R-K >= 2.7 mag^ and < J-K >= 0.9 
mag (Schombert et al. 1993; Silva & Elston 1994) and are in 
agreement with evolutionary synthesis models (e.g. Poggianti 
1997). 

Color profiles for the eight resolved galaxies are presented in 
Figure 4. We used R-band data from the HST survey (Scarpa 
et al. 2000a; Urry et al. 2000) and the present K-band results to 
compute these profiles. The BL Lac host galaxies show a clear 
gradient in their colors, being bluer toward the periphery. This 
trend was previously found by Kotilainen, Falomo, & Scarpa 
(1998b) and confirmed by Scarpa et al. (2000b) with additional 
data for more objects. We computed an average slope of the 
profiles to be A(R - K) /Alog r = -1 .27 ± 0.73 (Table 3). Com- 
pared with results from the previous NIR studies, this average 
color gradient appears, however, to be much steeper than ob- 
served in R-H. The magnitude of the average gradient values 
is reduced by only about one-half even when the three galax- 
ies exhibiting extremely large gradients hovering around -2 are 
not considered. The large gradients we computed may be an 
artifact due to the quite small effective radii we fit to our NIR 
data - our average K-band effective radius is 4.2 kpc, compared 
with ^10 kpc in the R-band (e.g., Falomo & Kotilainen 1999; 
Urry et al. 2000), although we can not be certain of the true 
cause. For comparison, normal non-active elhptical galaxies 
show much less pronounced gradients, but with the same sign, 
with typically A(y-K)/Alog r = -0.16 (Peletier, Valentijn, & 
Jameson 1990), and -0.26 (Schombert et al. 1993). 





5 




1 1 1 1 I < 1 1 1 1 1 
0158+002 






0331-363 : 


(0 


4 

3 










' I I w 

* ! i I i i i I J J i J 
















- K (rnag) 


5 
4 
3 


—+- 


H — I — I — — I — I — I — I — j — I — I — I — i 
0347-121 






— I — I — — I — I — I — I — — I — I — I — ^ 
0350-371 

"■ill, , 


K 


Z 












(mag) 


5 
4 


■—+- 


H — . — . — — . — . — . — . — 1 — I — I — . — i 
0521-365 






— I — I — — I — I — I — I — — I — I — I — ^ 
2005-489 : 




3 














2 












jnag) 


5 
4 


-—1- 


H — . — . — — . — . — 1 — I — 1 — I — I — 1 — h 
2143+070 






— 1 — 1 — 1 — — . — . — I — 1 — — I — I — \ — ^ 
2356-309 : 




3 














a 













10 20 10 2' 

Radius (kpc) Radius (kpc) 



Fig. 4. — R-K color profiles for the eight resolved host galaxies cal- 
culated using the R-band HST survey data (Urry et al. 2000; Scarpa et al. 
2000a) and the K-band data from the current work. The profiles indicate that 
the galaxies are bluer further from the nuclei. We find an average slope of 
A{R-K)/Alog r = -1.27 ± 0.73. See table 3 for individual values and sec- 
tion 3.2 for a discussion. 

al. 1993; Bressan, Chiosi, & Tantalo 1996). 
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3.3. Notes on the Host Galaxies of Individual Sources 

PKS 0138-097: The host galaxy of this high redshift (z = 
0.733) BL Lac object remains unresolved despite several at- 
tempts to separate the two components (e.g. Heidt et al. 1996; 
Scarpa et al. 1999a, 2000a). We masked out its three brightest 
companions before extracting the radial profile, but included 
the region in which the faintest companion (D) is embedded 
since it was found to contribute ^ 1% of the total flux at K'- 
band (Heidt et al. 1996). We found that our K-band profile 
is consistent with a 14.6 mag PSF, very near the K' value of 
14.4 mag found by Heidt et al. (1996). These values are nearly 
two magnitudes fainter than this object's historical K-band peak 
(Fan & Lin 1999). 

MS 0158+002: The host galaxy was clearly resolved by HST 
(Scarpa et al. 2000a; Urry et al. 2000) and ground-based obser- 
vatories (Wurtz, Stocke, & Yee 1996; Falomo & Kotilainen 
1999) as an elhptical with Mr = -23.0 mag, jig = 21.9 mag 
arcsec"^, and = 9.3 kpc (ffir values). There is a close com- 
panion to the N-W which we masked out in our analysis of the 
radial profile. We clearly detected the host (the first time in a 
NIR band), which appears to be a typical elliptical galaxy. 

MS 0331-363: An HST snapshot (Scarpa et al. 2000a; Urry 
et al. 2000) resolved the host as a giant (r^ = 18.7 kpc) elliptical 
galaxy with Mr = -24.3 mag and /i^ =22.1 mag arcsec"^. We 
have made the first detection of the host galaxy in the NIR and 
it is the largest galaxy (r^, = 8.5 kpc) of those we observed. 

lES 0347-121: Falomo et al. (2000) found a very circular 
elliptical galaxy hosting this BL Lac object from a 2-D struc- 
tural analysis of the HST survey data (Scarpa et al. 2000a; Urry 
et al. 2000). The //5rdata gave best fits of Mr = -23.2 mag and 
fie = 20.6 mag arcsec"^ which are consistent with ground-based 
values (Wurtz, Stocke, & Yee 1996; Falomo & Kotilainen 

1999) . However, the half-light radius obtained by HST (r^ = 5.3 
kpc) is about twice larger than the ground value (Wurtz, Stocke, 
& Yee 1996; Falomo & Kotilainen 1999). We resolve the host 
galaxy clearly, which is the first NIR detection of the host. Our 
value of r,, = 4.6 kpc is shghtly smaller than the HST value but 
is consistent with the observation by Scarpa et al. (2000b, their 
Figure 6) that BL Lac hosts tend to be bluer towards their outer 
regions. 

MS 0350-371: This BL Lac's host galaxy is a normal el- 
liptical as determined from 2-D analysis of the HST image by 
Falomo et al. (2000). The HST data are best fit by an ellipti- 
cal galaxy with Mr = -23.4 mag, /i^ = 20.8 mag arcsec"^, and 
re = 6.5 kpc (Scarpa et al. 2000a; Urry et al. 2000). We detected 
the host galaxy in two NIR bands - the first detections of the 
host at these wavelengths - and derived a NIR color {J—K = 0.7 
mag) very near the average in our sample. 

MS 0419+194: The host has been resolved in the optical 
band as an elhptical galaxy (Wurtz, Stocke, & Yee 1996; 
Falomo & Kotilainen 1999; Scarpa et al. 2000a; Urry et al. 

2000) . However, at z = 0.514, the host galaxy was unresolved 
in our mildly deep observation (27 min). Our analysis show 
that the 1-D K-band profile is consistent with a 16"' magnitude 
point source, very near the 15.4 mag observed by Gear (1993) 
during 1992. 

PKS 0521-365: This object contains a well-studied opti- 
cal/radio jet extending to the N-W (see Scarpa et al. 1999b, 
and references therein) which was recently detected in the X- 
rays (Birkinshaw, Worrall, & Hardcastle 2002). Its host was 
resolved in early ground-based imaging (Falomo 1994; Wurtz, 
Stocke, & Yee 1996), which are consistent with subsequent 



results (Mr = -23.3 mag, /i^ = 19.9 mag arcsec ^, and = 
4.1 kpc; Scarpa et al. 2000a; Urry et al. 2000). H-band imaging 
from the ground (Kotilainen, Falomo, & Scarpa 1998b) shows 
an intermediate galaxy absolute magnitude to our J and K re- 
sults (Mh = -25.4 mag), but for a larger galaxy (r^ = 3.8 kpc, 
versus our range of 1.5 to 2.5 kpc). Falomo et al. (2000) found 
this to be an average elhptical galaxy from a 2-D analysis of the 
HST data. 

H 1914-194: This is a little studied BL Lac object with 
no redshift determined. HST results show niR = 17.0 mag, 
/Xg = 24.7 mag arcsec"^, and r^ = 1" A for an elliptical host, but 
this was a rare instance in which a disk profile fits the //ST data 
equally wefl (Scarpa et al. 2000a; Urry et al. 2000). We fit a 
slightly smaller host (r^ = 4") in both of our K-band images, 
which are the first detections of the galaxy in the NIR. 

PKS 2005-489: The large round elhptical host of this weU- 
known blazar has been much studied. Our H-band results are 
remarkably consistent (within Icr) with those obtained by Koti- 
lainen, Falomo, & Scarpa (1998b). We resolved the bright 
galaxy in all three NIR bands even with short four to twelve 
minute exposures. Best fits from HST data (Mr = -23.9 mag, 
lie = 21.3 mag arcsec"^, and r^ = 10.5 kpc; Urry et al. 2000; 
Scarpa et al. 2000a; Falomo et al. 2000) show a bright elliptical 
host with derived parameters that agree with previous ground- 
based imaging (Stickel, Fried, & Kuhr 1993; Falomo 1996). 
However, the effective radii determined by the optical imag- 
ing is consistently larger than values derived in the NIR (Koti- 
lainen, Falomo, & Scarpa (1998b); this work) - up to three 
times larger than the J-band value, and five times larger at K. 
Our observations of this blazar were coordinated near a multi- 
wavelength monitoring campaign (Perlman et al. 1999; Taglia- 
ferri et al. 2001) in which both short and long-term X-ray out- 
bursts were observed. We found no significant variability of the 
nucleus between the nights in our NIR bands. 

MS 2143+070: The host galaxy of fliis BL Lac object has 
been well studied with HST and ground-based observatories. 
Scarpa et al. (2000b) derived niH = 15.3 mag from NICMOS 
data, which is consistent with the ground-based value derived 
by Kotilainen, Falomo, & Scarpa (1998b). These H-band val- 
ues are intermediate to our derived J and K galaxy magnitudes, 
as expected, and the effective radii we determined (6.0 and 7.1 
kpc, at K and J, respectively) are also within the range of H- 
band derived values from the ground (5.5 kpc) and space (7.6 
kpc). An I-band HST detection of the host by Jannuzi, Yanny, 
& Impey (1997) was confirmed (Urry et al. 1999) and com- 
bined with V-band data from the HST survey (Scarpa et al. 
2000a; Urry et al. 2000) and found to have aV-I color of 1.61 
mag. The HST survey results show Mr = -23.7 mag, He = 21.7 
mag arcsec"^, and r^ = 10.6 kpc (Scarpa et al. 2000a; Urry et 
al. 2000), assuming V -R = 0.86 mag, which agree well with 
ground-based results (Wurtz, Stocke, & Yee 1996; Falomo & 
Kotilainen 1999). 

PKS 2240-260: The redshift of 0.774 for this object is based 
on two weak features in the spectra and awaits confirmation 
(Stickel, Fried, & Kiihr 1993). The host was unresolved in the 
snapshot survey (Scarpa et al. 2000a; Urry et al. 2000) and 
remains unresolved by us, even with our fairly deep (42 min) 
image. This is the first pubhshed attempt to resolve the host in 
the NIR. We find a K-band magnitude for the BL Lac of 13.3 
mag that is near its historical maximum in this band (Fan & Lin 
1999). It has been seen to be as faint as 14.8 mag at this band 
(Gear 1993). 
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H 2356-309: The host is resolved as an elliptical in the HST 
data with Mr = -23 .2 mag, fie = 21.1 mag arcsec"^, and = 7. 1 
kpc (Scarpa et al. 2000a; Urry et al. 2000) in general agreement 
with previous ground-based results (Falomo 1991). Falomo et 
al. (2000) found that this is a normal elliptical from a 2-D anal- 
ysis of the HST data and also found a close 1".2 companion 
to the S-W. This companion is indistinguishable in our images 
and can contribute only < 1% to the total flux of the galaxy 
(Falomo et al. 2000) so was not masked. We imaged the galaxy 
in three NIR bands and our two epochs of H-band results agree 
well within our errors, and with results from NICMOS imag- 
ing {Mh " -25.5 mag and /ite = 18.5 mag arcsec"-^; Scarpa et al. 
2000b) except that they fit a larger galaxy (r^ = 6.5 kpc versus 
our derived value of 4.2 kpc) with smaller uncertainty. 

4. SUMMARY AND CONCLUSIONS 

We have successfully separated the NIR host galaxy emission 
from the bright nuclei of eight BL Lac objects. These objects 
are at known, and relatively low (z < 0.308) redshifts. Three 
targets at z ^ 0.5 were unresolved, and one marginally resolved 
BL Lac is at an unknown redshift. The hosts studied here are 
large elhptical galaxies with average luminosities comparable 
to brightest cluster member galaxies and over one magnitude 
brighter than normal ellipsoids (see section 3.1). This is consis- 
tent with findings established in previous H-band (Kotilainen, 
Falomo, & Scarpa 1998b; Scarpa et al. 2000b) and optical 
waveband studies (e.g. Falomo 1996; Wurtz, Stocke, & Yee 
1996; Scarpa et al. 2000a; Urry et al. 2000). Our additional data 
allowed us to better constrain a Kormendy relation for BL Lac 
objects in the NIR which is remarkably similar to that of normal 
ellipticals and in general agreement with observations of radio 
galaxies and other types of AGN (e.g. Taylor et al. 1996, see 
section 3.1). 



Our observations allowed us to investigate colors within the 
NIR-band for the first time in five galaxies, as well as pro- 
viding additional optical-to-NIR colors for eight objects. We 
find average R~K and J-K colors which are indistinguish- 
able from those measured in normal elliptical galaxies (sec- 
tion 3.2), typical of old stellar populations at the redshifts of 
the sources (Poggianti 1997; Scarpa et al. 2000b). The color 
profiles which utilized R-band HST observations show a bluer 
trend away from the galaxy centers, although the magnitude 
of this gradient is not definitive. Taken together, the presently 
available NIR data show that in all respects we could investi- 
gate, the BL Lac host galaxies are luminous ellipticals showing 
old stellar populations, and lends further support to a picture in 
which all elliptical galaxies can experience a phase of nuclear 
activity and that the galaxies do not know about the nucleus. 
Considering recent results on black hole demography (see Urry 
2003, for a recent summary), which shows that most probably 
all ellipticals harbor a super-massive black hole in their center 
(e.g. Gebhardt et al. 2000; Ferrarese et al. 2001), such a finding 
is not very surprising. 
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Table 1 



Journal of the Observations 



Name 


Date 


Filter 


Exposure* 


Sky Brightness 


FWHM 


Resolved? 




(in 1998) 




[minutes] 


[mag arcsec"^] 


[arcsec] 




PKS 0138-097 


Sep 29 


K 


54 


13.4 


0.79 


No 


MS 0158+002 


Sep 30 


K 


72 


13.4 


0.93 


Yes 


MS 0331-363 


Sep 29 


K 


36 


13.4 


0.83 


Yes 




Oct 1 


K 


18 


13.5 


1.05 


Yes 


lES 0347-121 


Sep 30 


K 


54 


13.3 


0.87 


Yes 


MS 0350-371 


Sep 30 


K 


18 


13.2 


0.93 


Yes 




Oct 1 


K, J 


36, 36 


13.5, 16.2 


0.92, 1.10 


Yes 


MS 0419+194 


Sep 28 


K 


27 


13.5 


0.93 


No 


PKS 0521-365 


Octl 


K, J 


28, 18 


13.3, 15.8 


1.01,0.91 


Yes 


H 1914-194 


Sep 28 


K 


63 


13.4 


0.95 


Yes 




Sep 30 


K 


54 


13.2 


0.72 


Yes 


PKS 2005-489 


Sep 29 


K, J, H 


4.5, 12, 4.5 


13.2, 15.8, 13.6 


0.78,0.89, 0.81 


Yes 




Octl 


K, J, H 


4, 12, 7.5 


13.2, 15.4, 13.6 


0.93,0.93,0.88 


Yes 


MS 2143+070 


Sep 29 


K 


54 


13.2 


0.89 


Yes 




Sep 30 


J 


54 


15.8 


0.92 


Yes 


PKS 2240-260 


Octl 


K 


42 


13.4 


0.93 


No 


H 2356-309 


Sep 28 


K, J, H 


54, 54, 27 


13.8, 16.7, 14.3 


0.82,0.88,0.86 


Yes 




Sep 29 


H 


54 


14.1 


0.86 


Yes 



Note. — 



*Total integration time for each night. 



Table 2 

BL Lac and Host Galaxy Near-Infrared Properties 



Name 


z 


Filter 




/T-con-C'' 


(c) 


(e) 
^PSF 


(c) 


/if) 


M 


PSF 




€ 








[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag/asec^] 


[arcsec] 


[mag] 


[mag] 


[kpc] 


PKS 0138-097 


0.733 


K 


0.02 




14.6 










-29.3 






MS 0158+002 


0.299 


K 


0.02 


-0.65 


14.9 


15.9 ±0.6 


14.6 ± 0.5 


17.0 ± 3.6 


0.8 ± 1.3 


-25.7 


-26.3 


4.7 ± 7.7 


MS 0331-363 


0.308 


K 


0.01 


-0.65 


14.7 


16.3 ± 0.8 


14.5 ± 0.3 


18.1 ±2.7 


1.4 ± 1.7 


-25.4 


-26.5 


8.5 ± 10.3 


lES 0347-121 


0.185 


K 


0.04 


-0.48 


14.5 


15.8 ± 0.5 


14.2 ± 0.3 


16.2 ± 1.1 


0.6 ± 0.3 


-24.7 


-25.8 


2.5 ± 1.3 


MS 0350-371 


0.165 


K 


0.01 


-0.44 


13.9 


15.4 ± 0.5 


13.7 ±0.2 


17.3 ± 1.5 


1.2 ±0.8 


-24.8 


-26.1 


4.6 ±3.1 






J 


0.03 


0.09 


15.3 


17.2 ± 1.1 


14.9 ± 0.4 


17.5 ± 2.6 


1.0 ± 1.2 


-23.0 


-25.4 


3.8 ±4.6 


MS 0419+194 


0.512 


K 


0.18 




16.0 










-27.1 






PKS 0521-365 


0.055 


K 


0.03 


-0.18 


11.3 


12.4 ± 0.3 


11.4 ±0.2 


14.7 ± 0.5 


1.0 ±0.2 


-25.3 


-26.1 


1.5 ±0.3 






J 


0.07 


0.03 


12.7 


14.1 ± 0.1 


12.7 ±0.1 


16.9 ±0.3 


1.7 ±0.2 


-23.6 


-25.1 


2.5 ± 0.3 


H 1914-194 




K 






12.7 


12.8 ± 0.1 


14.2 ± 0.4 




4.0 ± 1.9 








PKS 2005^89 


0.071 


K 


0.01 


-0.22 


10.2 


10.6 ± 0.2 


11.1 ±0.2 


14.4 ± 0.3 


1.0 ±0.1 


-27.6 


-26.9 


1.9 ±0.2 






H 


0.02 


0.00 


10.7 


11.1 ± 0.1 


11.7 ± 0.2 


15.7 ± 0.5 


1.5 ± 0.3 


-27.1 


-26.5 


2.8 ± 0.6 






J 


0.03 


0.04 


11.3 


11.7 ± 0.1 


12.6 ± 0.2 


17.1 ± 0.4 


2.0 ± 0.3 


-26.6 


-25.7 


3.7 ± 0.6 


MS 2143+070 


0.237 


K 


0.04 


-0.57 


14.4 


15.7 ± 0.5 


14.3 ±0.3 


17.7 ±2.6 


1.2 ± 1.4 


-25.4 


-26.2 


6.0 ±7.1 






J 


0.11 


0.11 


15.9 


16.8 ± 0.5 


15.8 ±0.5 


18.8 ±2.2 


1.4 ± 1.4 


-24.3 


-25.4 


7.1 ±7.1 


PKS 2240-260 


0.774te) 


K 


0.01 




13.3 










-30.8 






H 2356-309 


0.165 


K 


0.01 


-0.44 


13.8 


14.5 ± 0.2 


14.1 ±0.2 


17.5 ± 1.6 


1.1 ±0.8 


-25.7 


-25.6 


4.2 ±3.1 






H 


0.02 


0.01 


14.3 


15.1 ±0.3 


14.5 ±0.3 


17.4 ±2.2 


1.1 ± 1.1 


-25.1 


-25.7 


4.2 ±4.2 






J 


0.03 


0.09 


15.1 


15.9 ± 0.4 


15.2 ±0.4 


18.0 ± 1.6 


1.1 ± 1.7 


-24.3 


-25.1 


4.2 ± 6.5 


K-band median^''^ 


0.175 










15.6 


14.2 






-25.4 


-26.2 


4.4 



Note. — Interstellar extinction corrections using R-band extinction coefficients (Ag) from Urry et al. (2000) and appropriate NIR relations from 
Cardelli, Clayton, & Mathis (1989): Aj/Ar = 0311, Ah/Ar = 0.234, and Ak/Ar = 0.150. 

^'') K-corrections for first-ranked elliptical galaxies from Neugebauer et al. (1985) citing private communication from S.E. Persson (1984), interpolated to the 
redshifts of our targets. 

Observed total (Tot), BL Lac nucleus (PSF), and (Host) galaxy magnitudes. The host galaxy magnitudes were integrated to infinite radius so will appear 

in many cases to exceed the observed total magnitudes. 

Surface brightness at the half-light (effective) radius (re) with corrections for extinction, K-correction, and cosmological dimming applied. 
Bulge scale-length (effective radius), i.e., the radius at which half of the light in the galaxy is enclosed. 

Absolute magnitudes for the BL Lac nucleus (PSF) and Host. Both were extinction corrected but only the galaxy values were K-corrected. The BL Lac 
nuclei are assumed to have flat power law spectra, therefore have negligible K-corrections. 
Uncertain redshift. 
Statistics only for resolved objects. 
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Table 3 
Host Galaxy Colors 



Name 




j_K(") 


AiR-K)/Alogr '•'''> 




[mag] 


[mag] 


[mag] 


MS 0158+002 


3.3 




-2.00 


MS 0331-363 


2.2 




-1.10 


lES 0347-121 


2.6 




-1.06 


MS 0350-371 


2.7 


0.7 


-0.86 


PKS 0521-365 


2.8 


1.0 


-1.86 


PKS 2005^89 


3.0 


1.2 


-2.37 


MS 2143+070 


2.5 


0.8 


-0.33 


H 2356-309 


2.4 


0.5 


-0.57 


Median 


2.7 


0.8 


-1.08 



Note. — Integrated rest frame colors, extinction and 
K-corrected. R-band values are taken from Urry et al. 
(2000). 

Color gradient measured from the data presented in 
Figure 4 which used our K-band results and the R-band data 
from the HST survey (Scarpa et al. 2000a; Urry et al. 2000). 



